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1. Introduction 

The electromagnetic form factor of the pion G K (Q 2 ) is a simple and yet interesting quantity to 
explore hadronic structure using lattice QCD. At small momentum transfer, e.g., at Q 2 < 0.3 GeV 2 , 
the form factor is well parametrized in terms of the mean-square charge radius whose experimental 
value is (r 2 ) exp = 0.452(1 1) fm 2 . Since this is a well measured quantity, lattice calculation of (r 2 ^ 
is a good check for the validity of lattice QCD. 

In the framework of ChPT up to NLO, (r 2 ) was calculated by Gasser and Leutwyler in their 
famous work [|l|] : 

i / \ i m 2 \ 

<?)*j®#u> = -j. +w 2 + W 2lo8 W) ■ (L1) 

There is no m 2 suppression factor in front of the logarithm in this formula, and thus it predicts 
a rapid growth of the mean-square charge radius at small pion mass. There have been a number 
of efforts in the lattice QCD community to verify this behavior including recent studies with dy- 
namical quarks [§, ||, ||, ||] with the pion mass as light as about 300 MeV. However, the results 
for the mean-square charge radius at those pion mass range were significantly smaller than the 
experimental value, and in most cases these studies invoked the NNLO of ChPT to verify that the 
extrapolation toward the physical pion mass yields the charge radius consistant with experiment. 
Therefore a calculation closer to the physical point is needed to clarify the situation. 

Recently PACS-CS collaboration generated gauge ensembles of 2+1 dynamical flavors toward 
the physical point using the 0(a) improved Wilson-clover quarks and Iwasaki gauge action on a 
32 x 64 lattice volume with a lattice spacing of 0.0907(13) fm [Bj. Hopping parameters were 
chosen such that the pion mass covered the range from m n 702 MeV down to m n rj 156 MeV. 
It was shown that the spectrum quantities can be measured close to the physical point. This work 
encourages us to proceed to a project of calculation of the pion form factor on the PACS-CS gauge 
ensembles. This paper is a progress report of our work. 



2. Twisted boundary condition 

The minimum non-zero quark momentum for the periodic boundary condition on an 32 3 x 
64 lattice with a 2GeV inverse lattice spacing is about 0.4GeV. To probe the region of smaller 
momentum transfer one can apply the technique of twisted boundary condition [Q, ||, ^, [RJj. If one 
imposes the boundary condition given by 

y/(x + Lej) = e 2nid < yf{x) , j = 1,2,3 (2. 1) 

on quark fields where L denotes the spatial lattice size and ej the unit vector in the spatial j-th 
direction, the spatial momentum is quantized according to 

2Krij 2nd; 

Pi = —T + —r> 7 = 1,2,3. (2.2) 

In this way one can explore arbitrarily small momentum on the lattice by adjusting the value of 
twist dj. In practice we transfer the twist from the quark sector to the gluon sector by an U(l) 
transformation given by 

yr(x) — > U(0,x)yr(x) = e 2ni ^ e ^V(x). (2.3) 
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(a) Pion effective energies. 



(b) Continuum dispersion relation. 



Figure 1: Check of validity of the twisted boundary condition at K s = 0.1364, K ul j = 0.13700 where m n 
702 MeV with 40 configurations measured. 



Quark propagators are solved with the periodic boundary condition but on the gluon fields twisted 
by the U(l) transoformation above. 

To check that the term — j- 1 acts as a true physical momentum, we carried out a test on the 
relativity dispersion relation of the pion. Of the two valence quarks inside the pion, we twisted 
one quark with a twist angle 6 = (6,6,6) and left the other untwisted. In Figjlja) we plot the 
effective energy for several values of the twist angle and integer momenta at the hopping parameters 
K s = 0.1364, K lld = 0.13700 where m„ an 702 MeV. The propagator is fitted over t = 1 - 27 to 
extract the energy E(p). The error is estimated by the jackknife method with the bin size of 100 
trajectories. The results are plotted in Fig|l](b), together with the expected behavior, 

2 



E (p) = (am K ) + 



7.K _ 2n ■ . . 

— nH 6 
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(2-4) 



2nd, 



which demonstrates clearly that the term — j 2 - acts as a true physical momentum 



3. Electromagnetic form factor of the pion 

The electromagnetic form factor of the pion is defined by 

'7l+(p')\V 4 \7l + (p)) = (pA + PA)Gn(Q 2 ), 



(3.1) 



where p = (p,p4) is the 4-momentum of the incoming pion, p' = (p',p' 4 ) that of the outgoing pion, 
Q 2 = —(p' — p) 2 the 4-momentum transfer, and V4 = ^uy^u — ^dy^d the electromagnetic current. 
Since we employ the local (non-conserved) current for V4, the bare matrix element needs to be 
renormalized by requiring G K (0) = 1. 

One can extract the form factor from a suitable ratio of two- and three-point functions. In this 
work, we use the following ratio which has the advantage of simultaneously reducing fluctuations 



3 



The pion form factor from 2+1 dynamical flavors lattice QCD 



Oanh Hoang Nguyen 




°- 5 T 5 Tin fs st3 °- 5 6 crh 0^2 0^3 ~^ 0^5 0^ 0^7 as 

time Cr 



(a) Ratio R(t) at various momenta transfer. (b) Q 2 dependence of Gn(Q) 

Figure 2: Pion form factor at m n rj 702 MeV. (a) the ratio R(z) as a function of the time for the current. The 
sink is fixed at tf = 24, and the fitting range of G n (Q 2 ) is from 8 to 15. (b) Fitted values of G X (Q 2 ) and a 
line from the monopole fit. 




(a) Ratio R(t) at two momenta transfer. (b) Q 2 dependence of Gn(Q) 

Figure 3: The pion form factor at m K « 296 MeV. The sink time is fixed at tf — 24. The fitting range for 
G K (Q 2 ) is from 9 to 14. 



and renormalizing the current, 

C 3 ^(/,f /; p,0;T) C 



C 3 Py,f /; p,0;T)C 2 PV,T) 2E n {p>) larger and t f G^Q 2 ) r (n2 , 



(3.2) 

where C 2 ^(/?,t) = t)7T + (^,0)} and C 3 P'(p',t f ;p,0;i;) = (k+ (p' ,t f )%{x)\n + {p are 

the two- and three-point functions, respectively. 

Measurements were made for five values of the degenerate up-down hopping parameter k U( i = 
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M n (MeV) 


#conf measured 


e 


0.1364 0.13700 


702 


40 


0.18423, 0.28112 


0.13727 


570 


40 


0.18467, 0.28265 


0.13754 


411 


47 


0.18585, 0.28672 


0.13770 


296 


160 


0.18814, 0.29450 



Table 1: Simulation parameters. We adjust 9 values to have the same two smallest Q values at every m K 

{0.13700, 0.13727, 0.13754, 0.13770,0.13781}. The hopping parameter of strange quark is fixed 
at K s = 0.1364. The pion mass varies from M K ps 702 MeV down to M K w 156 MeV. We fix the 
final pion at zero momentum p' = and vary that of the initial pion p. The twist technique is 
applied to the quark running from source to the current. Two values were chosen for the twist 
angle 6 = (6,6,6) such that the smallest 4-momentum transfer of the current took the value 
<2 2 (GeV) = 0.01841 or 0.04237. Adding integer momenta, we then collected data for Q 2 in the 
range 0.1 GeV < Q 2 < 0.7 GeV. We employ an exponentially smeared source and local sink. The 



three-point function was calculated by the conventional source method [ |11| , [12 ]. The simulation 
parameters are listed in Table |T| 

In Fig^| we depict data for the pion form factor at pion mass m n 702 MeV. To extract the 
form factor, we fit the plateau of the ratio R(z) by a constant. The fitting range should be chosen 
around the symmetry point between the source and the sink, with additional considerations on the 
time interval required for the pion state to become dominant. In our case the sink is fixed at tf ■ = 24, 
thus the fitting should be symmetric around t = 12 if the sink and the source are of the same kind. 
However, since we use a smeared source and the point sink, we shift the fitting range one time unit 
to the source, and choose z = 8 to 15. 

Figure ^ shows similar data at the pion mass m n « 296 MeV. At current statistics, we have 
reasonable signal for only two smallest momenta transfer Q = {0.01841,0.04237}. 

4. Mononople ansatz and mean-square charge radius of the pion 

The experimental pion form factor follows the monopole form suggested by the vector domi- 
nance model 

Our data confirms this trend as illustrated in Fig.||. we then fit our data to the monopole form from 
which we extract the mean-square charge radius of the pion 

V) = 6 dQ 2 lff=0- ( 4 " 2 ) 

In Fig.|] we plot our initial results (filled circles) as a function of pion mass squared, together 
with the experiment value (burst) and those from latest researches: the only work for 2+1 dynamical 
flavor QCD is from RBC/UKQCD Collaboration at a single pion mass (open cirlce) using domain- 
wall quark acton, and others are for 2 dynamical flavors from ETMC Collaboration (twisted mass 
QCD; triangles) and from JLQCD (overlap action; inverted triangles). 
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5. Conclusions 



Figure 4: (r 2 ) dependence on M\ 



We have presented our initial results of a lattice calculation of the pion electromagnetic form 
factor in 2+1 dynamical flavor QCD with the O(a) improved Wilson-clover quarks and Iwasaki 
gauge action. Our current data of (r 2 ) at small pion masses agree with the recent data from other 
groups, and (V 2 ) show an increasing behavior toward the physical value as M n decreases. 

At present our data only goes down to M n rj 296 MeV whereas the PACS-CS gluon ensemble 
reaches M K m 156 MeV. Reducing the pion mass to 156 MeV, we have observed an increasingly 
larger fluctuation in the two- and three-point functions, and this trend worsens for larger momenta. 
For having better statistics, we are currently exploring various techniques. One of the possibilities 
is changing the set up to the channel where \p'\ = \p\. In this channel, the two-point function in 



the ratio ( p.2[ ) drops out, leaving just the ratio of three-point functions. This ratio of the three-point 
function with non-zero momentum transfer over that with zero momentum transfer seems to give 
signals with smaller statistical error compared with the ratio we used in the present calculation. We 
are also considering to use the random noise technique to reduce computer time for calculation of 
the two- and three-point functions. 
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